Spontaneous polarization and piezoelectricity in boron nitride nanotubes 
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Ab initio calculations of the spontaneous polarization and piezoelectric properties of boron nitride 
nanotubes show that they are excellent piezoelectric systems with response values larger than those of 
piezoelectric polymers. The intrinsic chiral symmetry of the nanotubes induces an exact cancellation 
of the total spontaneous polarization in ideal, isolated nanotubes of arbitrary indices. Breaking of 
this symmetry by inter-tube interaction or elastic deformations induces spontaneous polarization 
comparable to those of wurtzite semiconductors. 
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I. INTRODUCTION 

Piezo- and pyroelectric materials for modern techno- 
logical applications should display an excellent piezo- 
electric response, combined with high mechanical stabil- 
ity and low environmental impact. Existing materials, 
which can be broadly divided into the families of ce- 
ramics and polymers, can only partially fulfill the above 
requirements. Lead zirconate titanate (PZT) ceramics, 
for example, are strong piezo- and pyroelectricsi^ but, 
unfortunately, they are also brittle, heavy and toxic. 
On the other hand, polymers like polyvinylidene fluoride 
(PVDF) are lightweight, flexible and virtually inert, but 
their polar properties are an order of magnitude weaker 
than those of PZT. 3 In this paper, we examine sponta- 
neous polarization and piezoelectricity in boron nitride 
nanotubes (BNNTs) in order to estimate their poten- 
tial usefulness in various pyro- and piezoelectric device 
applications, and to understand the interplay between 
symmetry and polarization in nanotubular systems. 

BNNTs, broadly investigated since their initial predic- 
tion^ and succeeding discovery^ are already well known 
for their excellent mechanical properties. 6 However, un- 
like carbon nanotubes, most of BN structures are non- 
centrosymmetric and polar, which might suggest the ex- 
istence of non-zero spontaneous polarization fields. Re- 
cently, these properties have been partially explored by 
Mele and Krai, using a model electronic Hamiltonian. 7 
They predicted that BNNTs are piezo- and pyroelectric, 
with the direction of the spontaneous electric field that 
changes with the index of the tubes. The ab initio calcu- 
lations presented in this paper provide a much fuller de- 
scription and show that BNNT systems are indeed excel- 
lent lightweight piezoelectrics, with comparable or better 
piezoelectric responce and superior mechanical proper- 
ties than in piezoelectric polymers. However, contrary to 
the conclusions of Ref. 0, our combined Berry phase and 
Wannier function (WF) analysis demonstrates that elec- 
tronic polarization in BNNTs does not change its direc- 
tion but rather grows monotonically with the increasing 
diameter of the tube. Furthermore, the electronic and 



ionic spontaneous polarizations in BNNTs cancel exactly 
and these systems are pyroelectric only if their intrinsic 
helical symmetry is broken by, e.g., inter-tube interac- 
tions or elastic distortions. 

The rest of this paper is organized as follows: Sec. [H] 
briefly reviews the formulation of the modern polariza- 
tion theory in terms of Berry phases or Wannier func- 
tions. It also presents the details of the numerical tech- 
niques that were used to compute polarization. In Sec. lIIII 
we discuss the results and the complementary nature of 
the two techniques to compute the spontaneous polariza- 
tion. Finally, Sec. IIVI presents the summary and conclu- 



II. COMPUTATIONAL METHODS 

A. Modern theory of polarization 

The problem of computing polarization in materials is 
very subtle and is best approached by the "Berry-phase" 
method, introduced only a decade ago by Vanderbilt and 
King-Smith, 8 and Resta»2> Within this approach, the po- 
larization difference between two states of a system is 
computed as a geometrical quantum phase. In practice, 
this difference, AP = p( A i) - p( A o) ; can be obtained if 
one can find an adiabatic transformation A from one state 
to the other that leaves the system insulating. In the 
spirit of Ref. IToL p( A ) can be split into two parts: P} ^ 

and P^ , corresponding to the ionic and electronic con- 
tributions respectively. In the case of paired electron 
spins, the expression for the total polarization of the sys- 
tem can be written as follows: 
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where V is the volume of the unit cell, Z T and r T are the 
charge and position of the r-th atom in the cell, and uq- 
are the occupied cell-periodic Bloch states of the system. 
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For the electronic part, an electronic phase ip a (Berry 
phase) defined modulo 27r can be introduced as 

<pg>=VG a .P$>le, (2) 

where G a is the reciprocal lattice vector in the direction 
a. Similarly, one can construct an angular variable for 
the ionic part, called in what follows the "ionic" phase, 
so that the total geometrical phase is 

$W = ^zWG a .rW+,W. (3) 

r 

The total polarization in the direction a becomes 

= eQWRc/V, (4) 

where R a is the real-space lattice vector corresponding 
to G a , (R a ■ G a ) = 1. 

Alternatively, the electronic polarization of a system 
can be expressed in terms of the centers of charge of the 
Wannier functions of its occupied bands&S 

*s* = -|e/ 'i^ (A) wi a * = -|e^ a) >, (5) 

i i 

where the Wannier function (WF) W t (r) is constructed 

from the Bloch eigenstates of band i using the uni- 
tary transformation 

Wi(r) = / e ikr u ik (r)dk, (6) 

and (r^ A -* ) is the center of charge for the WF . How- 
ever, because of the arbitrariness in the choice of the 
phases of the Bloch orbitals (non-uniqucness of trans- 
formation JOJO) there is no unique representation of the 
WFs of a given group of bands. In our approach, we 
employ an algorithm that has been recently developed 
by Marzari and Vanderbilt^i which exploits the freedom 
in transformation © to construct WFs that are as lo- 
calized as possible. This is achieved by minimizing the 
sum of the quadratic spreads of the Wannier probability 
distributions | W}(r)| 2 , 

M 

^EKr-D-fo) 2 ], (7) 
i=i 

where the sum is over an isolated group of bands. The 
maximally localized WFs generated by this procedure are 
real, apart from an overall phase factor. 

In both methods presented above, -Pjj can be ob- 
tained only modulo 2eR/V due to the arbitrariness in 
the choice of the phases of the Bloch functions. How- 
ever, the difference in polarization AP is well defined 
if |AP e ;| <§C \2eR/V\. The same indetermination issues 
apply to ifil 



B. Calculations 

In computing the spontaneous polarization as the dif- 
ference between a polar BNNT and a nonpolar reference 
state, a natural choice for the nonpolar state is a nan- 
otube of the same geometry, but with boron and nitrogen 
atoms substituted by "pseudo carbon" atoms, which are 
50% boron and 50% nitrogen. The adiabatic transforma- 
tion is then defined by a "virtual crystal" procedure, in 
which parameter A corresponds to the content (in atomic 
%) of a site that is transformed from pure boron to the 
nonpolar reference state (vice- versa for nitrogen sites). 

We used an ab initio multigrid-based total-energy 
method, employing a real-space grid as a basis^ for all 
the Berry phase calculations presented here. The Ceper- 
ley- Alder— form, parametrized by Perdew and Zunger^4 
was used for the exchange-correlation energy functional 
in the local density approximation. The norm-conserving 
pseudopotentialsi^ for all the elements, including "vir- 
tual" ones, were generated by the f hi98PP packagsi^ uti- 
lizing the Klcinman-Bylandcr formulation— 

To isolate the contribution of individual nanotubes, we 
performed polarization calculations for periodic crystals 
of noninteracting (i.e., positioned sufficiently far apart) 
nanotubes in hexagonal and tetragonal arrangements. 
The electronic structure calculations were carried out us- 
ing two special fc-points along the T-A direction in the 
hexagonal or T-Z direction in the tetragonal Brillouin 
zones. The fc-space integration to compute ifz was done 
on a string of 20 fc-points uniformly distributed along the 
same direction and shifted to avoid the T-point. The in- 
ternal consistency of our approach was checked against 
the results obtained using the ABINIT code^ for a few 
selected systems with excellent agreement. 

Because of the different alignments of the polar bond 
with respect to the nanotube axis, we anticipate that the 
symmetry of the nanotube will play an important role in 
determining the magnitude of the spontaneous polariza- 
tion field. In particular, since the zigzag geometry maxi- 
mizes the axial dipole moment, we expect to observe the 
strongest effects in (ra, 0) nanotubes. 



III. RESULTS AND DISCUSSION 

A. Berry-phase method 

The ionic part of the polarization in zigzag BNNTs, 
presented in Fig. ^ is large and directly proportional to 
the nanotube's index. This is in contrast, for instance, to 
the corresponding wurtzite III-V and II- VI systemsjiSiSl 
where the spontaneous polarization can be viewed as the 
difference between the polarizations of the wurtzite (po- 
lar) and zincblende (nonpolar) geometries. Since these 
configurations become geometrically distinct only in the 
second shell of neighbors, their ionic phases are very 
close. The major contribution to the spontaneous polar- 
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ization in wurtzite materials is then due to the difference 
between the electronic polarizations (which are 0.04-0.08 
C/m 2 ), while in BNNTs both the ionic and the electronic 
contributions are essential. 
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FIG. 1: Ionic- phase difference between the polar and non- 
polar configurations for zigzag nanotubes; the ionic phase of 
the nonpolar configuration is set to zero. Inset: ionic phases 
wrapped into the [—tt, tt] interval. Phases are given in units 
of n. 

The ionic phase differences A.tpi on between the polar 
and nonpolar configurations of zigzag nanotubes were 
evaluated via the virtual crystal approximation. The 
inset in Fig. shows the results obtained by a simple 
lattice summation over the ionic charges (the first term 
in Eq.|2J), with the phases translated into the [—tt, tt] in- 
terval. The phases plotted in the main graph were "un- 
folded" by eliminating all the 2tt discontinuities and set- 
ting the phase of the nonpolar reference configuration to 
zero. For the unfolded phases, as the diameter of a nan- 
otube increases, i.e., as another hexagon is added around 
the circumference of the tube, the ionic phase goes up by 
7r/3, so that the total ionic phase for a (n, 0) BN nan- 
otube amounts to rnr/3. 

In Fig.[3we show the electronic-phase differences Aip e i 
between the polar and nonpolar configurations for zigzag 
nanotubes. These data suggest a natural division of the 
nanotubes into three families with different Atp e f. tt/3 
for n = 31 — 1, — 7r/3 for n = 31 + 1, and — tt for n = 31, 
where I is an integer (2i which is similar to the result 
obtained by Mele and Krali However, the existence of 
such three classes of behavior is surprising, given that 
the ionic character of the electronic charge density (as- 
sociated with the B-N bond) does not change with the 
nanotube index. Additionally, there is an important dif- 
ference between our results and those of Ref. 0, where 
the electronic polarization of heteropolar nanotubes was 
studied within a simple 7r-orbital tight-binding (7T-TB) 
approximation. In Ref. |3, the "n = 3?" family has a zero 
electronic phase instead of —tt. 
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FIG. 2: Electronic-phase differences between the polar and 
nonpolar configurations for zigzag nanotubes. 



This discrepancy is due to the ambiguity of the 
definition of electronic polarization as a multivalued 
quantity^ which can assume a lattice of values corre- 
sponding to Berry phases that differ by arbitrary multi- 
ples of 2tt. Unlike the ionic phase model, where discon- 
tinuities in <^i on (A) can be easily monitored, Berry phase 
calculations always produce phases that are smoothly 
folded into the [—tt, tt] interval and cannot be extrapo- 
lated. To obtain an unambiguous determination of the 
spontaneous polarization of BNNTs of arbitrary diame- 
ters, one has to compute the polarization in a different 
way, using the centers of charge of the WFs of the occu- 
pied bands (Eq. Note that this approach does not 
solve the problem of branch indetermination, since while 
Berry phases are defined modulo 2tt, Wannier centers are 
defined modulo a lattice vector R. However, by shifting 
the indetermination from the phase to the lattice vector, 
we are able to map the electronic polarization problem 
onto a simple electrostatic model, where the unfolding of 
the electronic phase is straightforward^ 



B. Maximally localized Wannier functions 

The results of the maximally localized WF calculations 
for BNNTs are summarized in Fig. El where examples of 
the WFs for C and BN zigzag nanotubes of arbitrary 
diameter are shown, together with a schematic drawing 
that illustrates the shift of the Wannier centers in the 
adiabatic transformation from C to BN. Since 



(BN- 



= _| E(r (-) 



(8) 



the magnitude of the shift of the centers is directly pro- 
portional to the electronic polarization of the BNNT with 
respect to the nonpolar CNT. 

The cr-band WFs are centered in the middle of the 
C-C bonds in carbon nanotubes, while they are shifted 



4 




FIG. 3: Upper panel: Examples of Wannier functions (WF) 
of the a and tt occupied bands of C (left panel) and BN (right 
panel) nanotubes. Lower panel: Schematic positions of the 
centers of the Wannier functions in C and BN hexagons, and 
the projections of the tt WF onto the nanotube axes. The 
positions of the centers of a WF are indicated by circles, and 
those of 7r by squares. The direction of the shifts of a and 7r 
WF in an adiabatic transformation from C to BN is indicated 
by arrows. The projections of shifts of the a WF cancel, so 
that the a WF do not contribute to polarization (see text). 



towards the cations in BN nanotubes because of the dif- 
ferent electronegativities of B and N atoms. Since these 
shifts have the same magnitude along each of the three 
bond directions, the vector sum of all shifts is zero (see 
bottom panel of Fig. [3J , and the a orbitals do not con- 
tribute to the total polarization of the system. The tt- 
band WFs are centered on the cations in BNNTs, while in 
CNTs they have a peculiar V-shape, with centers some- 
what outside of the C-C bond. The sum of the shifts of 
the 7r-band Wannier centers is non-zero only for the axial 
component, which means that the electronic polarization 
in BNNTs is purely axial. 

The bottom panel of Fig.[3]shows the projection of the 
tt WF centers onto the axis of the tube. The projections 
of the centers have an effective periodicity of half of the 
axial lattice constant c, which leads to the indetermina- 
tion of the electronic phase by multiples of tt. Moreover, 
the WF description allows for an unambiguous unfolding 
of the electronic phase. In analogy to the ionic phase, 
we find that each individual hexagon carries a phase of 
— 7r/3, leading to a total electronic phase of — n7r/3 for a 
(n, 0) nanotube. This result demonstrates that the direc- 
tion of the electronic polarization in a BNNT is specified 
by the orientation of the B-N bond and does not oscil- 
late in direction with the nanotube diameter, contrary 
to the model Hamiltonian predictions^ We should point 



out that the Wannier function results are completely con- 
sistent with the Berry-phase calculations, since an elec- 
tronic phase of —mr/3 for any n can be folded, modulo 
tt, into the three families found previously. 

When we combine the results for the ionic and elec- 
tronic phases into a general formula for the phase of an 
arbitrary (n, m) BNNT, 



A$^ ot (n,m) = A<pl on (n,m) 



■ A(pf(n,m) 
n — m 



(9) 



we find that the two contributions cancel exactly and 
that the total spontaneous polarization in any BNNT 
is zero, i.e. the Wannier centers are arranged in such a 
way as to completely compensate the polarization due to 
ions. We have verified this result by two-point (A = 50 
and 100%) calculations of the Berry phase difference for 
a number of chiral nanotubes ((3,1), (3,2), (4,1), (4,2), 
(5,2) and (8,2)) and found an exact cancellation in all 
BNNTs, except for those narrower than approximately 
4 A, where a residual polarization is present as an effect 
of the very high curvature. In such nanotubes Wannier 
centers cannot fully compensate the ionic polarization, 
due to the severe distortion of the atomic bonds, which 
makes these systems weakly pyroelectric. For example, 
P = 0.11 C/m 2 in (3,1), 0.008 C/m 2 in (7,0) and 0.002 
C/m 2 in (12,0) nanotubes. 

The exact cancellation is a result of the overall chi- 
ral symmetry of the nanotubes which, although not cen- 
trosymmetric, are intrinsically nonpolar. Nevertheless, 
cancellation of ionic and electronic polarizations is exact 
only in the limit of an isolated BNNT. The spontaneous 
polarization in a nanotube bundle, where the chiral sym- 
metry is effectively broken, is different from zero. For ex- 
ample, in (7,0) bundles at equilibrium distance of 3.2 A 
P w 0.01 C/m 2 . However, in this case it is hard to 
estimate the separate contributions to polarization due 
to bundling, extreme curvature and elastic deformation. 
Although smaller than in polymers or PZT, this polariza- 
tion is comparable to some wurtzite pyroelectrics: e.g., 
P = 0.06 C/m 2 in w-ZnO^. 



C. Piezoelectricity 

The Berry-phase method can also be employed to com- 
pute piezoelectric properties of BNNTs, which are di- 
rectly related to polarization differences between strained 
and unstrained tubes. In the linear regime, the change 
in polarization due to strain can be decomposed into a 
sum of two terms: a uniform axial strain and a relative 
displacement of the two sublattices. It is therefore natu- 
ral to describe the geometry of a BNNT of a given radius 
in terms of an axial lattice constant c and an internal 
parameter u, where uc is the length of the vector con- 
necting the anion with the cation. With this choice, the 



5 



TABLE I: Piezoelectric properties of zigzag BNNT bundles!** 
The corresponding values for a few piezoelectric materials are 
listed for comparison. 



(n, m) 


diameter (A) 


Z* (e) 


|e 33 | (C/m 2 ) 


Ref. 


(5,0) 


3.91 


2.739 


0.389 




(6,0) 


4.69 


2.696 


0.332 




(7,0) 


5.47 


2.655 


0.293 




(8,0) 


6.24 


2.639 


0.263 




(9,0) 


7.04 


2.634 


0.239 




(10,0) 


7.83 


2.626 


0.224 




(11,0) 


8.57 


2.614 


0.211 




(12,0) 


9.38 


2.609 


0.198 




(13,0) 


10.16 


2.605 


0.186 




tu-AIN 




2.653 


1.50 


[19] 


w-ZnO 




2.11 


0.89 


[20J 


PbTiOa 






3.23 


\2] 


P(VDF/TrFE) 






« 0.12 


\3] 



axial piezoelectric polarization is 

dP 3 dP 3 

5P 3 = £33^3 = ~^(c- C ) + -q—(u- U ), (10) 

where the strain is 63 = (c — co)/co, and cq and uo are 
the equilibrium values of c and u. The only surviving 
piezoelectric strain tensor component is 

(°) _l ec o AT 7 * du m\ 
e 33 = e 33 ' + — NZ — , (11) 

where N is the number of B-N pairs in the supercell. 
Here, 

efi = Ct) dP 3 /dc (12) 

is the "clamped-ion" piezoelectric constant (representing 
the effect of strain on the electronic structure), and 

Z* = {V/eNc„)dP 3 /du (13) 

is the axial component of the Born dynamical charge ten- 
sor. Both polarization derivatives were computed as fi- 
nite differences, changing c or u by ±1%. The parameter 
£ = codu/dc, describing the change in the bond lengths 
under axial strain, was obtained by rescaling c together 
with the associated components of ionic coordinates, and 



then relaxing the geometry of the system. For all the sys- 
tems considered below, the value of £ is approximately 
the same and equal to -0.085. 

We have calculated the piezoelectric properties for var- 
ious bundles comprised of zigzag BNNTs with individual 
diameters ranging from 3.9 to 10.2 A. These results are 
summarized in Tableland compared to a few well-known 
piezo- and pyroelectric materials. While the piezoelectric 
constants of zigzag BNNTs are modest when compared 
with inorganic compounds, they are still substantially 
larger than those in the PVDF polymer family. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the spontaneous po- 
larization and piezoelectric properties of BN nanotubes 
using state-of-the-art ab initio methods. Our calcula- 
tions demonstrate the complementary nature of Berry 
phase and Wannier function analysis, and show that a 
real-space description is necessary to unravel the Berry 
phases in complicated cases. The results suggest that 
BNNTs are excellent nonpolar piezoelectrics that exhibit 
substantially higher strain response than polar polymers. 
Moreover, we have shown that, contrary to the previous 
expectations, ideal non-interacting nanotubes are effec- 
tively nonpolar due to their intrinsic chiral symmetry, 
which leads to a total cancellation between the ionic and 
electronic polarizations. Breakage of this symmetry, as 
in the simple case of interacting nanotubes in a bundle, 
induces spontaneous polarization fields that are compa- 
rable to those of wurtzite semiconductors. Due to their 
piezo- and pyroelectric properties, BNNTs are excellent 
candidates for various nano-electro-mechanical applica- 
tions. 
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